Linear temperature dependence of conductivity in Si two-dimensional electrons near 

the apparent metal-to-insulator transition 
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In a high mobility two-dimensional electron system in Si, near the critical density, n c = 
0.32 x 10 n cm -2 , of the apparent metal-to- insulator transition, the conductivity displays a lin- 
ear temperature (T) dependence around the Fermi temperature. When ao, the extrapolated T = 
conductivity from the linear T-dependence, is plotted as a function of density, two regimes with 
different ao(n) relations are seen, suggestive of two different phases. Interestingly, a sharp transi- 
tion between these two regimes coincides with n c , and ao of the transition is ~ e 2 /h, the quantum 
conductance, per square. Toward T = 0, the data deviate from linear <r(T) relation and we discuss 
the possible percolation type of transition in our Si sample. 

PACS numbers: 72.20.-i, 71.30.+h 



The nature of the ground state of an interacting two- 
dimensional electron system (2DES) in the presence of 
disorder is a long-standing problem in condensed mat- 
ter physics. More than a decade ago, an apparent 2D 
metal-to-insulator transition (MIT) was first reported as 
the density of a Si MOSFET is reduced through a char- 
acteristic density tu \m . Despite much research effort 
later in the field |2j, |3j, there remain several unsettled 
fundamental questions. For example, metallic behavior 
is observed for density n > n c . Does this metallic-like 
state persist down to T = and thus represent a true 2D 
metal? Besides, the 2DES eventually becomes insulating 
as ever-decreasing densities. Is this phenomenon related 
to a novel phase transition due to strong electron-electron 
(e-e) interactions or is it a mundane crossover just 
due to a complex combination of many well-understood 
physical mechanisms 

In an earlier publication, we reported the observation 
of the metallic behavior and the apparent 2D MIT in 
a high mobility Si quantum well |fj. The high-density 
metallic-like state and its response to an in-plane mag- 
netic field are emphasized in that paper. In this commu- 
nication, we focus on the transport properties for densi- 
ties close to the transition. Here, the conductivity dis- 
plays a linear temperature dependence near T = Tp (the 
Fermi temperature) and the slope is the same for all dif- 
ferent densities around n c . When the extrapolated T = 
conductivity, ao, of this T-dependence is plotted as a 
function of density, two regimes with linear ao vs n re- 
lations are readily seen. Interestingly, these two <7o(n) 
lines cross each other exactly at n c and ao at the cross- 
ing coincides with the quantum conductance, e 2 /h, per 
square. At low temperatures, T << Tp, the measured 
a{T) deviates from the linear T-dependence. We discuss 
the low-T behavior of our data within the percolation 
model. 

The experiments were performed on the 2DES 



in an n-type Si quantum well confined in a 
Sio.75Geo.25/Si/Sio.75Geo. 25 heterostructure. The 
2D electron density is tuned continuously by applying 
a front gate voltage to our field-effect transistor device. 
Details of the growth and the sample structure can 
be found in Ref. |6|. Standard low-frequency (~ 
7Hz) lock-in techniques were used to measure the 2D 
resistivity p. At T ~ 300mK and zero gate voltage, the 
2DES has a density n = 1.45 x 10 u cm -2 and mobility 
fi = 190, 000cm 2 /Vs. 

In Fig. 1(a), we reproduce selected T-dependence data 
p(T) from our previous paper Toward the T = 
limit, the apparent 2D MIT is clearly seen at the critical 
density n c = 0.32 x 10 u cm -2 , where dp/dT ~ for T < 
IK. The insulating behavior, dp/dT < 0, is observed for 
densities below n c and metallic-like behavior, dp/dT > 0, 
above n c . For further insights to the 2D MIT, we now 
focus on the data around n c and plot the inverse resis- 
tivity, or conductivity a, as a function of temperature in 
the density range of 0.27 x 10 11 < n < 0.38 x 10 n cm _2 in 
Fig. 1(b). In the low-T limit, the MIT is again observed 
in that da/dT changes sign as n changes through n c . At 
high temperatures when T approaches the Fermi temper- 
ature Tp, marked as short lines for each density, however, 
all a(T) curves show roughly a linear T-dependence and 
bends slightly downward for T sufficiently higher than 

To illustrate this linear a(T) relation around T = Tp, 
we show the linear fits, a(T) = ao+"/T, for several densi- 
ties, where 7 is the slope and ao the linear extrapolation 
to T = 0. In the inset of Fig. 2, 7 is plotted as a function 
of the density. It is nearly constant, ~ 0.43 ± 0.01e 2 //i 
per Kelvin, for n < 0.38 x 10 11 cm -2 and then rapidly 
decreases for higher densities. We plot ao vs n in Fig. 
2 and two regimes with different linear ao{n) relations 
are readily identified. On the low-density insulating side, 
ao increases as increasing density at a rate of 11 e 2 /h 
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FIG. 1: (a) Selected data of p(T) in a high quality Si quan- 
tum well reproduced from Fig. 1 in Ref. [5]. (b) u(T) in 
the vicinity of the critical density. The electron densities, 
from bottom to top, are 0.268, 0.282, 0.296, 0.303, 0.318, 
0.322, 0.324, 0.328, 0.336, 0.344, 0.352, 0.358, 0.365, 0.374, 
and 0.380, in units of 10 11 cm -2 , respectively. The short, ver- 
tical lines mark the Fermi energy at each density. The lines 
on some curves are the linear fits around T — Tf- 



per 10 11 cm -2 . The slope roughly doubles to 26 e 2 /h 
per lO^cm -2 on the high-density side. Strikingly, the 
do (n) data show a sharp bend at the crossing of the two 
straight lines at n = 0.32 x 10 1:L cni _2 , which coincides 
with n c of the 2D MIT. And <tq at this density exactly 
equals the quantum conductance, e 2 /h, per square. For 
both regimes, we extrapolate the linear a® (n) to zero con- 
ductivity at m = 0.23 x 10 11 cm -2 on the low-density side 
and ri2 = 0.28 x 10 n cm -2 on the high-density side. 

Linear cr(T) relation has been observed before in p- 
type GaAs samples 0. There, 7 = dcr/dT - 3 to 5 
e 2 /h per Kelvin, compared to the much smaller number, 
~ 0.43e 2 //i per Kelvin, in our Si quantum well sample. 
This one order of magnitude difference might be related 
to the different Fermi temperatures in these two systems. 
As will be shown below, der/dT ~ 1/Tp. In our Si sample, 
Tp (~ 2 to 2.5K) is about a factor of 10 larger than that 
(~ 0.1 to 0.4 K) in the p-GaAs samples. As a result, 
der/dT in Si is expected to be ~10 times smaller than in 
p-GaAs. 

Theoretically, the linear cr(T) in the high temperature 
regime has also been addressed by several models. Das 
Sarma and Hwang [9J calculated the temperature de- 
pendence of a within a classical model of the screened 
charged impurity scattering. For T >> Tp, they ob- 
served that cr(T) ~ T/Tp oc T, consistent with our ex- 
perimental observation. Besides the screening model, the 
micro-emulsion model proposed by Spivak and Kivelson 
can also explain the linear (T) at high temperatures. 
In this model, near n c , the ground state of the 2DES con- 
sists of an electron liquid with Wigner crystal inclusions 
on the metallic side, and Wigner solid with liquid in- 
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FIG. 2: The T = conductivity, ero, extrapolated from the 
high-T linear fit in Fig. 1(b), is plotted as a function of the 
2DES density. The solid lines are linear fits on both low-n 
and high-n regimes and they extrapolate to <jo = at ni = 
0.23 x 10 n cnr 2 , and n 2 = 0.28 x lO^cm -2 , respectively. 
The dotted lines show that the two linear fits cross at n c and 
cro = e 2 fh. The inset shows the slope, 7, of the linear c(T) in 
Fig. 1(b), as a function of density, is roughly constant below 
0.38 x 10 11 cm -2 , and decreases rapidly for higher densities. 



elusions on the insulating side. At high temperatures, 
Wigner crystal droplets melt and the system behaves 
classically. In analogy to the physics of 3 He near the 
crystallization pressure lQf , the viscosity of the electron 
liquid, which is directly proportional to the resistivity, is 
inversely proportional to T, resulting in a linear cr(T), 
again consistent with our observation [IH . 

Quantitatively, however, neither of the two above mod- 
els can explain the experimental observation that dcr / dT 
is nearly independent of the electron density around 
n c . Under the temperature dependent screening model, 
dcr/dT is proportional to 1/Tp or n _1 As a result, 
the slope should decrease by a factor of ~ 1.5 when the 
2DES density increases from 0.27 to 0.38 x 10 n cm" 2 . On 
the other hand, the micro-emulsion model shows dcr/dT 
~ n 2 |lflj . Consequently, in the same density range, the 
slope should change roughly by a factor of 2. In con- 
trast, dcr/dT is nearly constant and ~ 0.43 ± 0.01e 2 //i 
per Kelvin in our measurements. So far, it is not known 
what is responsible for this inconsistency between our 
experimental result and the theoretical predictions. 

Having discussed the linear T-dependence of the data, 
we need to address some puzzling aspects of the overall 
results. First, the examined specimen has very high elec- 
tron mobility and the n c of the 2D MIT is by far the 
lowest among all the Si-based samples. The e-e interac- 
tion parameter r s at the transition density is ~ 10, i.e., 
the Coulomb energy E c exceeds the Fermi energy Ep by a 
factor of 20, after taking into account the two-fold valley 
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FIG. 3: Percolation plot of <r(0.3K) vs n/n p — 1 in a log-log 
scale. Three trial densities at 0.23, 0.28, and 0.32x 10 n cm~ 2 , 
indicated by different symbols, are used as n p in the plot. 
For n p = 0.28 x 10 11 cm -2 , excluding the two lowest density 
points, a linear fit (solid line) is obtained over two decades in 
a. And the slope corresponds to an exponent of a = 1.31. 

degeneracy in (001) Si 2DES. Consequently, even though 
the system behaves classically for T ~ Tp, the 2DES is 
still strongly correlated since E c is the dominant energy 
scale here. Second, if the high-T physics underlying the 
linear a(T) behavior persists to the T = limit, our 
<Jo(n) data show two distinct regimes, possibly suggest- 
ing two different electronic phases below and above n c . 
The transition between the two density regimes is sharp, 
and occurs almost exactly at n c . a® at this transition 
point is very close to e 2 /h, the unit of quantum conduc- 
tance, per square, manifesting a possible quantum nature 
of this transition. The reduction of the rate, dao/dn, of 
the density dependence, from ~ 26e 2 /h per 10 11 cm~ 2 
above n c to ~ lie 2 /h per 10 11 cm -2 below n c , clearly 
indicates that in the low-density regime the 2D electrons 
are less likely to become localized by reducing the 2DES 
density. Considering the large r s at n c , one might specu- 
late that strong e-e interactions help to prevent electrons 
from being localized at low densities. Indeed, in a re- 
cent publication 0], Shi and Xie showed that, with e-e 
interactions taken into account, the 2DES becomes less 
localized compared to a non-interacting system at the 
same density. However, the appearance of two density 
regimes is still unexplained because a smooth evolution 
should be expected from their calculations. 

We now turn to the low-T limit of the measured <j(T) 
in Fig. 1(b). Instead of following the trend from high 
temperatures, the conductivity deviates from the linear 
T-dependence toward the T — ground state configura- 
tion. It has long been suggested that the apparent MIT 
observed here might be of a density inhomogeneity effect, 
and belongs to a general class of 2D percolation problem 



[l3lll^ |. In this picture, at very low densities, the 2DES is 
macroscopically inhomogeneous and first forms isolated 
puddles by occupying the low potential "valleys" . As n 
increases, the area of the electron puddles increases and 
at the percolation threshold density n p , some puddles 
are connected, giving rise to a conducting path through- 
out the sample. Experimentally, this percolation-type 2D 
MIT is supported by various measurements, e.g., scan- 
ning near-field optical microscopy (SNOM) [Iff, scanning 
single-electron transistor (SET) microscopy bjji trans- 
port 0], and surface-acoustic wave (SAW) ^| experi- 
ments. 

According to the percolation model, the conductivity 
of the 2DES follows the scaling function a ~ (n/n p — 1 ) Q , 
and the exponent a in a classical percolation transition is 
4/3 [14]. In Fig. 3, we show the measured a at T = 0.3K 
as a function of n/n p — 1 in a log-log scale. Since n p 
is not known beforehand from finite temperature mea- 
surements, three trial densities of n p are used, in unit 
of 10 11 cm -2 , n\ = 0.23 (extrapolation to <tq = from 
the low-n regime), n 2 — 0.28 (extrapolation to do = 
from the high-n regime), and n c = 0.32 (critical density 
of the observed 2D MIT). For n p = n 2 = 0.28 x 10 11 
cm~ 2 , except for some deviation at very low densities, a 
good linear fit with a = 1.31, close to 4/3 in the classi- 
cal model, can be obtained over two decades in a. It is 
not known, at this stage, whether this coincidence is ac- 
cidental or the two densities are actually deeply related. 
On the other hand, for n p — ri\ or n c , the fittings to a 
power law behavior are poor. We nevertheless empha- 
size two important points when applying the percolation 
model to our data. First, the power law of conductiv- 
ity in the percolation model only holds for densities very 
close to n p , or n/n p — 1 << 1 [13]. Second, percolation 
transition is essentially a zero-temperature phase transi- 
tion. As seen in Fig. 1(b), cr{T) does not saturate at our 
lowest measured T = 0.3K. In this regard, it is neces- 
sary that further measurements be carried out at lower 
temperatures. 

In summary, in a high quality Si quantum well speci- 
men, near the apparent 2D metal-to-insulator transition, 
a linear temperature dependence of conductivity is ob- 
served at T around Tp on both sides of n c . When <7o, 
the extrapolation of this linear a(T) to T = 0, is plot- 
ted as a function of density, two regimes with different 
do vs n relations are readily seen. Interestingly, the two 
linear ao(n) regions cross almost exactly at n c , and Co 
at the crossing point is e 2 /h, the quantum conductance, 
per square. We also show that the measured a(T) at our 
low-T limit can be fitted by a percolation scaling func- 
tion a ~ (njn p — 1) Q when n p — n%, the extrapolation of 
ao(n) to T = on the high-density side. 
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